Quarter-filled extended Hubbard model with alternating transfer integral: 
Two-dimensional Ising transition in the ground state 
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We study the one-dimensional quarter-filled extended Hubbard model with an alternating transfer 
integral. In the strong-dimerization limit the charge part is described by the quantum Ising model 
which shows the two-dimensional Ising criticality at the self-dual point, and it is naturally connected 
to the double- frequency sine-Gordon theory in the weak dimerization. Treating low-lying excitations 
in finite-size systems, we numerically determine a phase boundary between two types of 4fcp density- 
wave states and clarify the ground-state phase diagram. Further, we refer to its relevances to the 
charge-ordered phase observed in the charge-transfer organic salts. 
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The organic materials described by the chemical for- 
mula (TMTSF) 2 X (X=PF 6 , C10 4 , etc) and (TMTTF) 2 X 
(X=PF6, AsFg, etc) form a class of the quasi one- 
dimensional (ID) conductors; a large number of investi- 
gations on these materials have been accumulated in the 
literature^ While the various types of electronic phases, 
e.g., spin/charge-density- wave (SDW/CDW), the spin- 
Peierls and the superconducting states have been ob- 
served in the low-temperature region, newly discovered 
charge-ordered (CO) phase in (TMTTF) 2 X exhibiting an 
anomaly in the low-frequency dielectric constant^ and 
the charge disproportionation in the NMR studies^ has 
received intensive current interest. Although the stabi- 
lizations of these phases at finite temperature resort to 
interchain couplings, it is believed that the intrachain 
interaction effects play a leading role to describe them. 

For the study of the CO phase, the ID quarter-filled ex- 
tended Hubbard model (EHM) with an alternating trans- 
fer integral has been employedi 4 i 5 i 6 i 7 i 8 iV° H = H 1 + H 2 
with 

Hi = E -* t 1 - ( c Ls+m + H - c -) . (!) 

H 2 =Y^ ( Un 3,T n J,l + VrijUj+x) , (2) 
3 

where Cj s annihilates an s-spin electron (s ="f or J,) on 
the jth site and satisfies the periodic boundary condi- 
tion c L+1 s = c 1 s (j € L is an even number). 

The number operators are defined as rij }S — cj s Cj S and 
n j = n jA^ n j,i- The parameters U and V taking positive 
values stand for the onsite and nearest-neighbor Coulomb 
repulsion. The dimerization parameter 6 shows the alter- 
nation in the transfer integral of the molecular chains (we 
set t = 1 in the following). For the theoretical descrip- 
tions of the ID electrons, the Tomonaga-Luttinger liq- 
uid (TLL) picture has been widely adoptediii Since TLL 
consists of the massless charge and spin parts both con- 
trolled by the Gaussian fixed point [the conformal field 
theory (CFT) with the central charge c = 1], it is impor- 



tant to understand its instabilities. In particular, the CO 
transition may be related to the crossover of the critical- 
ity embedded in the renormalization group (RG) flow, 8 
which is one of the typical instability of the c =1 CFT. 

In this paper, we present the numerical calculation 
results of the ground-state phase diagrams of H. Our 
method being deeply connected to the instability will 
be explained briefly. Further, in the strong-dimerization 
limit (5 = 1), we show that the charge part can be de- 
scribed by the so-called quantum Ising chain, which is 
complementary to the bosonization argument and gives 
us an exact limiting condition of the phase boundary line. 
Since the region with sufficiently large Coulomb repul- 
sions is relevant to the CO transition, an occurrence of 
the phase separation or the transition to the supercon- 
ducting phases is outside of our research scope. 

Let us start with the description of the low-energy 
physics in the weak-coupling region, where the bosoniza- 
tion method provides a reliable approach, i.e., linearizing 
the dispersion at Fermi points ±fcp = ±nn/2a (electron 
density n = N/L = 1/2) and applying the method, we 
can obtain an effective Hamiltonian. For the present case, 
according to the recent research results ^LLiS we can use 
the following expression: H —>H = TL p + Ha with 
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-g p smV8(f> p + g 1 / 4 cos2V8(j) p ) (3) 



where the operator 9 P is the dual field of <p p satisfying the 
commutation relation [<f> p (x), d y 6 p {y) / it] = iS(x — y) and 
parameters K p and v p are the Gaussian coupling and the 
velocity of the charge excitation, respectively^ A ben- 
efit to use the bosonized expression is now clear, i.e., 
since the spin-charge separation occurs in TL and Tier is 
the SU(2) critical Gaussian model in the present case^^ 
we can concentrate on the charge part H. p which takes 
a form of the so-called double-frequency sine-Gordon 
(DSG) model. In uniform case (5 = 0), the 8fcp-Umklapp 
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scattering with g 1/4 cx U 2 {U - 4V) (Refs. Il2ll4ll5ll6j) 

brings about the Berezinskii-Kosterlitz-Thouless (BKT) 
transition, and then the charge part becomes massive 
for large values of the Coulomb interactions. For the 
BKT transition point, values in the strong coupling limit 
are known as V*(U — ► oo) = 2 and U*(V — ► oo) = 
Further the estimations for the intermediate re- 
gion are availableiiSiiSiiS In the case of nonzero dimer- 
ization (5 ^ 0), the scaling dimension of the "half- filled 
Umklapp scattering" term with g p cx US[1 — A(U — 2V)] 
(A is a constant )LA& on the Gaussian fixed point is 
small (.T4B = 2K p ) enough to bring about the second- 
order phase transition for V < V*(U) 1 which is accom- 
panied by the divergent correlation length of the form 
\ cx < 5- 1 /(2-2ifp)j - 2 a For y > V *(U), since the charge 
gap may survive in a weak-dimerization region, the tran- 
sition point S P (U, V) takes nonzero values depending on 
U and V, and more importantly, the universality of the 
transition is changed. Recently, Tsuchiizu and Orignac, 8 
on the basis of the DSG theoryj2I argued that the charge 
part on S p (U,V) [V > V*(U)} is renormalized to the 
2D-Ising fixed point with c = ^ (i.e., the fixed point 
with lower symmetry), which is in accord with Zamolod- 
chikov's c-theorem 2 ^ (see also Refs. |23| and |24|) . Then, 
the critical line corresponds to the phase boundary and 
satisfies a condition S P (U, V — > V*(U)) \ in the weak- 
dimerization region. To characterize the phases, we shall 
use the CDW and the bond-order- wave (BOW) order pa- 
rameters with the 4fcp wave vectort 1 ^ 

04C cx cos V8<fi p , 04B oc sin VS(f> p . (4) 

Here, note that the expectation value of the 4/cp-BOW 
order parameter is finite, (04b) and (C>4c) = in the 
upper region of the boundary, but both of these are finite 
in the lower region (5^0). While this "mixed" state is 
basically the 4fcp-CDW phase, we shall use the double 
quotation marks "4/cp-CDW" to express this situation. 24 
On the other hand, another condition of the boundary 
can be found in the strong-dimerization limit (5 = 1). 
To derive an effective Hamiltonian, it is convenient to 
work with the orbital operators defined by d m ± s = 
( C 2m-M ± c 2m,s) I where d m l s annihilates an s-spin 
electron in the Z-orbital (/ = ±) on the mth unit cell 
(to e [l,L/2]). In this limit, Hi consists of a sum of 
the intracell electron hopping, which is diagonalized by 

using the operators as Hi = E m ,i s - 2l 4n i,a d m,l,s- For 
sufficiently large U and V, since the one-electron states 
\h s )m = dl nls \0) have a principal role to describe the 
mth unit cell in the quarter-filled ground state, and the 
Hamiltonian does not change the electron number in each 
cell, we shall introduce the pseudospin operators, 

Tm = \ d ln,l, S Mu' d rn,l>,s, (5) 

l,l',a 

acting on the orbital space as, for instance, T^|±, s) m = 
±i|±,s) m [t = (t 1 ,t 2 ,t 3 ); t % is the Pauli matrix]. Us- 



ing these, Hi — — 4T^. For H 2 , since the intra- 
cell Coulomb interactions are absent and the intercell 
Coulomb repulsion only remains in the restricted Hilbcrt 
space spanned by the direct product of one-particle states 
{cE>m|^s)m}, a straightforward calculation brings about 
the expression H 2 = J2 m ( — ^m^m+i + const). Now, 
since the Hamiltonian acts only on the orbital space, its 
eigenstate takes a form of the direct product of vectors in 
the spin and the orbital spaces as |$) = |spin) ® | orbital). 
Thus, assuming a certain spin configuration belonging to 
the 2 L / 2 -dimcnsional space for spins and restricting our- 
selves to the orbital (or charge) part, we see that the 
Hamiltonian H with S = 1 is reduced to the quantum 
Ising chain 2 ^ 

H P .S=1 = ^ (~ rT m ~ JT m T m+l) ( 6 ) 
m 

(r = 4, J = V). Note that this possibility was mentioned 
qualitatively in Ref.^J Then, the ground state of Eq. © 
is known to show the 2D-Ising criticality at its self-dual 
point T = J/2 (V — 8), which separates ordered ((T^) ^ 
0) and disordered ((T^) = 0) phases. The ordered state 
is realized via the breaking of the Z 2 symmetry (r 1 — > 
— r 1 ), and it is doubly degenerated, e.g., 

I ±rl ) =n^( d m,+ ± 4,-)|0) =II4m-l(cL)|0)(7) 

m m 

(we dropped the spin index). This expresses the 4fcp- 
CDW state with the perfect microscopic polarization, 
(ir 1 ^! ± r 1 ) = (±T%{n 2m -i - n 2m )\ ± r 1 ) = ±i 
On one hand, a disordered state is supported by the ex- 
ternal field in redirection, and an ideal one is given by 

I + - 114+1°) = II 71 (4m-l + 4n)lo>. ( 8 ) 

which expresses the 4fep-BOW state as expected. Here 
it is worthy of noticing that these states can be distin- 
guished by the expectation value of the twist operator 2 - 

3 

This quantity takes values z p = 1 for Eq. Q and 
z p = -[cos(27r/L)] L/2 for Eq. ©, so the sign of z p 
characterizes these two density- wave states (see below). 
Consequently, in the strong-dimerization limit, the or- 
bital degrees of freedom show the 2D-Ising type transi- 
tion between the "4fc F -CDW" and the 4fc F -BOW phases 
at V = 8, where U is irrelevant. Since this pseudospin 
representation is naturally connected to the bosonization 
picture in the weak couplings, 8 the phase boundary be- 
longs to the 2D-Ising universality and satisfies the limit- 
ing condition S P (U, V — > 8) /* 1, which provides a solid 
guide to investigations in the strong-dimerization region. 

Here, note that the qualitative estimation of the phase 
boundary might be possible in the weak- and strong- 
dimerization region^ To evaluate the entire phase dia- 
gram precisely, however, a numerical treatment of the 
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ID electron model is required. For this issue, recently the 
present authors have numerically treated the same insta- 
bility observed in the quantum-spin chain and interacting 
electron systems^ Therefore, we shall emplo y th e same 
approach to the present system (see also Ref. 28). Since 
there are two critical fixed points connected by the RG 
flow, a relationship between lower-energy excitations on 
these fixed points — the ultraviolet-infrared (UV-IR) op- 
erator correspondence — has essential significance in the 
investigationsA2i To see this, let us rescale phase fields 
and the Gaussian coupling as 2<fi p — > cf>, p /2 — > 9, and 
AK p — ► K ~ 1, which makes it possible to directly adopt 
our previous research^ With respect to c/>, the nonlinear 
potential density is given as — g p sin ^/2<j> + g\u cos \/S(f), 
and the order parameters as O4C oc cos y/2c/> (x4c = K/2) 
and O4B oc sin\/20 (^4B = K/2). Along the RG flow 
these operators on the Gaussian fixed point (UV) are 
transmuted to those on the 2D-Ising fixed point (IR) as 

4C ^M, Om^I + e, (10) 

where fj, is the disorder field (Z2 odd), and e is the en- 
ergy density operator (Z 2 even) with scaling dimensions 
= g and x e = 1, respectively. Since the dimeriza- 
tion 5 couples with O4B in the Hamiltonian J3J), a de- 
viation from the transition point S — 6 p (U,V) plays a 
role of the "thermal scaling variable" and brings about 
£ oc [6- 5 p {U,V)]- u with \jv = 2 - x t = 1. On one 
hand, the operator fi corresponding to O4C provides a 
most divergent fluctuation. 

Now, we shall explain our numerical procedure to de- 
termine the transition point. We shall focus our atten- 
tion on the level AE in finite-size systems which cor- 
responds to the operator 4C (taking the ground-state 
energy as zero). According to the finite-size-scaling ar- 
gument based on CFT, AE ~ 2irx4c/L on the UV fixed 
point)2S we can numerically obtain the level by using dis- 
crete symmetries of the lattice Hamiltonian in the di- 
agonalization calculations. Various excitations observed 
in TLL are characterized by a set of quantum num- 
bers for symmetry operations. With respect to 04c, it 
can be found in the subspace of the total spin S5, = 
and the space inversion P = — 1 (the boundary condi- 
tion is the same as that for the ground state) >i& Sup- 
pose that AE(U, V, S, L) is a level corresponding to O4C 
in the L-site system. Then, we numerically solve the 
phenomenological renormalization-group (PRG) equa- 
tion (L + 2)AE(U, V,S,L + 2) = LAE(U, V, S, L) with 
respect to S for given values of U and V, where the 
gap behaves as AE oc \/L [i.e., an L-dependent transi- 
tion point S P (U, V, L + 1) (see Fig. After evaluating 
S P (U, V,L + 1), we extrapolate them to the limit L — > 00 
using the formula S p (U, V, L) = 5 p (U, V) +aL~ 3 ^S where 
S P (U, V) and a are determined by the least-square-fitting 
condition. 

From the data of L = 12-20, we obtain the phase 
boundary S p (U, V) as shown in Fig.^ We can check that, 
for all values of U used here, the phase boundary lines 
converge to the point (V,S) = (8,1) with the 2D-Ising 





FIG. 1: The ground-state phase diagram of the quarter-filled 
EHM with an alternating transfer integral. The correspon- 
dence between marks and system sizes is given in the figure. 
The double circles show the limiting values, i.e., (V*(U), 0) at 
which criticality changes from the Gaussian to the 2D-Ising 
type, and (8, 1) the self-dual point. 
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FIG. 2: The 5 dependence of z p (L) at U = 16 and V = 4. 
The vertical dotted line indicates the transition point (i.e., 
the PRG result) 5 P (16, 4) ~ 0.12. Inset plots L dependences 
of S' p (U, V, L) (crosses) and the PRG data 5 P (U, V, L) (circles) 
with the fitting line. 

criticality. On the other hand, while the finite-size correc- 
tions to the boundary may be large in weak-dimerization 
region, the boundaries also show convergences to the 
BKT-transition points (V* (U),0). Next we demonstrate 
the S dependence of z p (L) in Fig. With the increase 
of S, z p (L) decreases and becomes negative [we denote 
the zero point of z p {L) as S'(U,V,L)]; this corresponds 
to the change of the center of mass as demonstrated in 
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the abovei2& However, unlike, for instance, the Gaussian 
transition, z p (oo) can take a finite value on the Ising tran- 
sition point, so S'JU, V, L) may not give an estimation of 
the transition point. In fact, the inset of Fig. [21 exhibits 
that S' p (U, V, L) may be extrapolated to a value different 
from the PRG result. On the other hand, Fig. also 
shows that there is a point S ~ 0.12 at which z p (L) is 
almost independent of L. This crossing point is expected 
to be a good estimator for the Ising transition point be- 
cause this is quite close to the PRG result even for small 
L. However, this issue remains as a future problem. 

Lastly, we shall refer to some implications of our study 
to the real materials. Besides the quantum-chemistry 
calculations the numerical estimations of the model 
parameters have been performed based upon the exper- 
imental dataA£ For example, the realistic values of the 
dimerization parameter and the onsite Coulomb repul- 
sion of (TMTTF) 2 PF 6 have been estimated as t 2 /t 1 ~ 
0.7, U/h ~ 7.0 (t lj2 = 1 ± S), but the value of V is 
still controversial (an uncertainty exists also in the value 
of t 2 /t^2£). Our numerical estimation of the transi- 
tion point using these values is V c /ti ~ 4.0, while gen- 
erally the mean-field-type calculations tend to predict 
somewhat smaller values due to an overestimation of V 
effectsAifi On the other hand, several values have been 
reported for this material, e.g., V/t\ ~ 2.8 (1.4) in Ref. 
0(Ref.@>, which is much smaller than the critical value, 



and thus predicts a uniform charge distribution (this con- 
clusion may not be changed even in smaller dimerization 
cases). However, (TMTTF) 2 PF 6 has the CO phase in the 
region above the lower-temperature spin-Peierls phase, 
and further it was theoretically suggested that a huge 
anomaly in the dielectric constant may reflect a nature 
of systems in the critical region^ This discrepancy may 
be attributed to many other interaction effects not in- 
cluded in the Hamiltonian. However, we think that since 
experimental findings seem to support the spin-charge 
separation with respect to the CO transition?^ the ID 
electron models are to provide a primary description of 
real materials. 

To summarize, we investigated the ground-state phase 
diagram of the ID quarter-filled extended Hubbard 
model with alternating transfer integral. Especially, the 
criticality on the phase boundary and the implication to 
the CO transition observed in the charge-transfer organic 
salts were mainly argued. 
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